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by 

Doreen  M.  C.  Walker 

SUMMARY 

The  rocket  of  Cosmos  268,  I969-20B,  entered  orbit  on  5  March  1969,  with 
an  initial  perigee  height  of  230  km  and  inclination  of  48.40°.  Accurate 
orbits  were  computed  at  RAE  from  all  available  observations.  Using  the  values 
of  perigee  height  from  the  RAE  orbit  and  decay  rates  from  Spacetrack  bulletins 
103  values  of  density  have  been  calculated  between  July  1969  and  February  1970 
On  three  occasions  when  geomagnetic  activity  was  strong  there  were  sudden 
increases  iri  density.  When  the  density  was  corrected  to  a  fixed  height,  the 
semi-annual  variation  was  apparent.  There  was  a  strong  minimum  in  .July  1969, 
a  maximum  in  October-November  1969  and  a  weak  minimum  in  January  1970. 
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INTRODUCTION 


The  air  density  at  heights  near  180  km  has  been  determined  from  the. 

orbit  of  1967-3! A  *,  between  July  1968  and  September  1969,  and  at  heights  near 

2 

150  km,  between  December  1969  and  August  1970,  from  1969-108A  .  These  two 

studies  leave  a  gap  in  time  between  September  and  December  1969  and  the 
present  study  was  undertaken  to  evaluate  the  variations  in  density  between  the 
two,  and  to  overlap  at  each  end.  Unfortunately  there  was  no  satellite  avail¬ 
able  with  a  perigee  height  below  200  km.  The  satellite  selected  was  the 
rocket  of  Cosmos  268,  1969-20B.  The  density  was  evaluated  over  the  last  seven 
months  of  t'ie  satellite’s  life,  from  July  1969  to  February  1970.  During  this 
time  the  perigee  height  decreased  from  218  km  to  155  km. 

Cosmos  268,  I969-20A,  waj  launched  on  5  March  1969,  and  its  rocket, 

19'69-20B,  entered  an  orbit  with  an  initial  perigee  height  of  230  km,  an  apogee 

height  of  2170  km,  and  inclination  48.40°.  This  rocket,  about  8  m  long  and 
3 

1.65  m  m  diameter  ,  was  a  good  object  for  optical  observation  from  Britain, 
and  because  of  its  interest  for  upper-atmosphere  research,  it  was  given  high 
priority  for  observing  by  t'nc  British  optical  and  radar  tracking  stations, 
including  the  Hewitt  camera  at  Malvern.  Accurate  orbits  have  been  computed 

4 

at  RAE  from  all  available  observations  ,  the  perigee  heights  having  an  average 
sd  of  about  200  m.  In  the  present  paper  these  perigee  heights  were  used, 
together  with  the  orbital  decay  rates  calculated  from  USAF  Spacetrack  elements: 
103  values  of  air  density  at  heights  between  185  km  end  253  km  were  obtained 
at  dates  between  17  July  1969  and  8  February  1970. 

2  ORBITAL  DECAY  RATE 

The  orbital  decay  rate  was  found  with  adequate  accuracy  from  the  USAF 

Spacetrack  5-card  elements,  of  which  114  sets  were  issued  during  the  seven 

1  2 

months.  As  in  previous  studies  ’  the  values  of  the  mean  anomalistic  motion 
n  (rev/day)  were  tabulated  and  differenced  to  give  values  of  An,  each  of 
which  was  divided  by  the  corresponding  time  interval  At  (day)  to  give 
values  of  n  *  An/At  (rev/ day  ).  Any  erroneous  value  of  n  causes  an  up- 
down  pattern  in  h,  i.e.  one  value  is  too  high  and  the  next  too  low,  or 
viae  versa.  After  inspecting  the  values  of  n  up  to  MJD  40608,  six  values  of 
n  were  eliminated  which  led  either  to  up-down  patterns  or  to  time  intervals 
At  less  than  24  hours;  one  more  value  of  n  was  omitted  over  the  final 
phase,  MJD  40608-MJD  40628,  because  of  an  up-down  tendency,  but  the  minimum 
limit  on  At  was  reduced  to  6  hours. 
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The  remaining  106  values  of  n  are  plotted  as  circles  in  Fig.l,  and  the 
values  from  the  RAE  orbits  are  plotted  as  triangles.  The  RAE  vaJues  a**  for 
specific  dates  while  the  n  values  from  Spacetrack  are  averages  over  a  time 
interval  At,  With  this  difference  taken  into  account,  the  two  sets  show 
excellent  agreement. 

3  PERIGEE  HEIGHT 


3. 1  The  ’perigee  height  parameter'  Q 

The  perigee  distance  a(l  -  e),  where  a  is  the  semi  major  axis  and  e 
the  eccentricity,  suffers  a  regular  oscillation  caused  by  the  effect  of  the 
odd  zonal  harmonics  in  the  Earth's  gravitational  potential.  For  1969-20B  at 
an  inclination  i  of  48.4°  this  perturbation"*  amounts  to  5.10  sin  km,  where 
u  is  the  argument  of  perigea.  In  the  Spacetrack  orbital  elements  this 
perturbation  is  taken  to  be  £.8  sin  i  sin  m  *  5.09  sin  w  km,  and  is  subtracted 
from  the  values  of  e.  Since  the  Spacetrack  perigee  distances  are  not  likely 
to  have  an  accuracy  better  than  0.5  km,  the  difference  of  0.01  km  is  not 
significant. 


It  is  useful  to  define  a  'perigee  height  parameter’,  Q  say,  which  is 
free  of  the  effects  of  the  odd  2onal  harmonics.  The  obvious  definition  to 
use  is 


ia(l  -  e)  -  R] 


Spacetrack 


or 


Q  »  [a(l  -  e)  -  +  5.10  sin  u 


(1) 


where  R  is  the  Earth’s  equatorial  radius  (6378.16  km).  The  quantities  given 
by  these  equations  should  be  identical,  apart  from  observational  errors,  and 
Q  should  decrease  slowly  under  the  combined  influence  of  air  drag  and 
luni-solar  perturbations. 

Fig. 2  gives  values  of  Q  from  the  Spacetrack  elements  (circles)  and  the 
RAE  orbits  (triangles).  The  RAE  values  all  have  standard  deviations  less  than  , 
350  m,  and  should  be  more  accurate  than  the  Spacetrack  values.  A  smooth  curve 
has  been  drawn  through  the  values  of  Q  in  Fig. 2,  giving  more  weight  to  the 
accurate  RAE  values,  but  the  real  variation  of  Q  with  time  will  not  be  smooth 
because  of  the  variations  in  n  apparent  in  Fig.l. 
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3.2  Calculation  of  perigee  height  y 


The  actual  perigee  height  y^  above  the  Earth  is  found  from  Q  by 
restoring  the  odd-zonal-harmonic  oscillation  (5.10  sin  c),  subtracting  the 


local  Earth  radius  at  latitude 


V 


which  differs  from  R  by 


2  2 

21.38  sin  £  =11.95  sin  m  kra,  and  adding  the  small  amount  by  which  the 

P  g  ,2 

actual  perigee  distance  differs  from  a(l  -  a),  namely  0.83  -  1.74  sin  . 

Thus  we  find 


km. 


y 


p 


Q  +  0.83  +  10.21  sin2  -  -  5.10  sin  - 
Q  +  q 


(2) 


The  variation  of  q  with  is  given  in  Fig.  3,  which  show's  that  the  combined 
effect  of  F,arLh  oblateness  and  odd  harmonics  can  alter  perigee  height  by  up  to 
16  km;  as  usual  the  maximum  perigee  height  occurs  at  w  -  270°. 


Fig. 4  show's  the  values  of  y^  for  1969-20B  as  given  by  equation  (2), 
with  Q  having  the  values  given  by  the  curve  in  Fig. 2. 

4  CALCULATION  OF  AIR  DENSITY 

The  air  density  at  a  height  $H*  above  perigee  was  found  from  the  usual 
equation 


0.157 ( lOOOn) 
106n2a 


-  8e  + 


7H* 

16ae 


^  +  £  si**2  1 


COS  2o;  I 


...  (3) 


where  H*  is  the  estimated  value  of  H  at  perigee  height,  and  may  be  in 
error  by  up  to  25%  without  incurring  an  error  of  more  than  1.2%  in  p^.  The 
assumed  variation  of  H  with  height  is  given  in  Fig. 5,  the  values  having  been 

7 

taken  from  the  latest  available  model  atmosphere  of  Jacchia  .  Equation  (3) 
is  valid  if  3H*/a  <  e  <  0.2,  conditions  that  are  satisfied  for  all  but  the 

last  three  values  of  n.  The  values  of  o.  in  equation  (3)  are  obtained  in 

3  A  2 

kg/m  if  H*  and  a  are  in  km,  the  area/mass  parameter  6  is  in  m  /kg, 

2 

•  n  is  in  rev/day  and  n  in  rev/day  . 

The  value  of  the  area/mass  parameter  6  is  not  known  accurately  for 

2 

-  Cosmos  316,  and  a  constant  value,  0,0246  m  /kg,  was  used,  corresponding  to  a 

2 

mass/area  ratio  of  83  kg/m  .  This  value  was  chosen  so  that  the  values  of 
density  at  the  average  of  the  semi-annual  variation  at  heights  between  180  and 
190  km  were  in  agreement  with  Jacchia’s  model2  for  the  appropriate  value  of 
exospheric  temperature,  1100  K, 


Equation  (3)  gives  values  of  oA»  the  air  density  at  a  height 
yA  *  y^  +  $H*.  In  examining  day-to-day  variations  in  density  we  need  to  con¬ 
vert  the  values  of  density  to  a  fixed  height,  y^,  say,  which  may  most 
appropriately  be  taken  as  the  mean  value  of  yA  over  a  suitable  time  interval. 
The  density  at  height  y.  is  calculated  from  the  equation 

15  iS 


=  exp 


(¥*) 


where  Hg  is  the  value  of  H  at  height  yB,'*  If  Hg  is  in  error  by  N%,  the 
use  of  equation  (4)  will  lead  to  an  error  of ■ approximately  [ (yA  -  yg)N/Hg|% 
in  c,,.  Since  H  may  be  in  error  by  up  to' perhaps  102  as  a  result  of 

D  D 

day-to-day  variations  in  H,  the  value  of  (y,  -  yn)/Hn  needs  to  be  kept  as 
small  as  possible.  It  was  therefore  necessary  to  split  the  total  time  interval 
into  several  parts.  The  breaks  were  made  at  MJD  40540,  40616  and  40624,  and 
the  values  of  the  various  parameters  during  the  three  time  intervals  are  given 
in  Table  1 . 


Table  1 

"he  four  sets  of  values  of  y  ,  H*,  yg  and  Hg 


Time  interval 
MJD 

Mean 

km 

y  H* 

km 

yB 

km 

hb 

km 

40419-40540 

225 

40 

245 

43 

40540-40616 

203 

34 

220 

38.5 

40616-40624 

185 

30 

200 

33.5 

40624-40626 

172 

26 

185 

29.5 

_ — . 

5  RESULTS 

5.1  Variation  of  density  with  time 


The  values  of  pg  obtained  are  plotted  as  points  joined  by  lines  in 
Fig. 6.  From  MJD  40419  to  MJD  40540  the  values  of  pg  are  for  a  height  of 
245  km,  i.e.  they  are  values  of  ^245*  *rom  40540  to  MJD  40616  the 
density  is  for  a  height  220  km,  but  the  quantity  plotted  is  0.538p22Q»  the 
numerical  factor  0.538  being  introduced  to  avoid  discontinuity;  from 
MJD  40636  to  40624  and  from  MJD  40624  to  40626,  the  values  plotted  are 
0.305p20Q  and  0.187pja5  respectively,  the  numerical  factor  being  included  for 
the  same  reason. 
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The  quantity  plotted  at  the  top  of  Fig. 6  is  the  sun-perigee  angle,  with 
the  local  time  at  perigee  marked  on  the  curve. 

5.2  The  day-to-night  variation 

In  Fig. 6  the  semi-annual  variation  in  air  density  is  somewhat  obscured 
by  the  day-to-night  variation,  and  a  correction  for  day-to-night  variation  must 
therefore  be  made. 

This  correction  was  calculated  by  using  Ref, 7.  From  Table  6  of  Ref. 7, 

values  of  p/p  are  plotted  against  T/Tq  for  various  heights  in  Fig, 7,  where 

p  is  the  value  of  density  for  a  given  temperature  T,  p  is  the  density  at 

temperature  1000  K  and  Tq  is  the  minimum  night-time  temperature,  taken  as 

870  K.  The  value  of  T/T^  appropriate  to  the  local  time,  solar  declination 

and  perigee  latitude  at  each  of  the  103  calculated  values  of  density,  was  then 

obtained  from  Table  1  of  Ref. 7;  and  hence,  knowing  the  perigee  heights  from 

Fig. 4,  values  of  p/p  were  found  from  Fig. 7.  These  are  plotted  against  MJD 

in  Fig. 8.  The  correction  for  the  day-to-night  variation  can  then  be  made 

by  dividing  each  calculated  value  of  density  by  the  appropriate  value  of 

p/p.  The  corrected  values  of  density  are  plotted  in  Fig. 9,  numerical  factors 

again  being  included  to  avoid  discontinuity.  The  quantities  plotted  at  the 

top  of  Fig. 9  are  the  solar  radiation  energy  on  a  wavelength  of  10.7  cm, 

as  measured  by  NRC,  Ottawa  and  the  daily  geomagnetic  index  A  as  giveh  by 

P  / 

the  Institut  fUr  Geophysik,  Gottingen,  plotted  with  a  12-hour  time  lag. 

5.3  Comparison  with  geomagnetic  and  solar  indices 

Between  July  1969  and  February  1970  there  were  three  occasions  when 
geomagnetic  disturbances  wera  strong  enough  to  give  values  of  greater 
than  40.  These  three  disturbances,  at  MJD  40430,  40495  and  40535,  are 
accompanied  by  marked  increases  in  air  density,  as  Fig, 9  shows. 

1  2 

Previous  studies  *  have  suggested  that  at  heights  below  200  km  there  is 
little  resemblance  be^reen  the  variations  in  air  density  and  in  the  solar 
10.7-cm  radiation,  v,r.-  h  is  a  rough  measure  of  the  extreme  ultraviolet  radia¬ 
tion,  The  height  of  I969-20B  is  slightly  over  200  km  until  the  last  15  days 
of  the  satellite's  life,  but  the  correlation  between  the  density  values  and 
SJ0  j  is  not  very  strong.  However  there  are  some  occasions  when  maxima  and 
particularly  minima  of  the  solar  activity  coincide  with  maxima  or  minima  in 
density.  For  example  the  minima  in  density,  from  Fig. 9,  at  MJD  40422,  40450, 
40478,  40510,  40565,  40595  and  40622  all  occur  when  the  values  of  SJQ  ?  are 
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at  minima.  There  are  also  three  increases  in  density  which  correspond  with 
rises  in  solar  activity,  namely  at  MJD  40460,  405*5  and  40600.  Of  these  three 
the  largest  is  that  around  MJD  4CI0,  amounting  to  about  40%  between 
MJD  40596  and  MJD  40603. 

5.4  The  semi-annual  variation 

If  the  peaks  of  density  at  the  times  of  geomagnetic  disturbance  are 
ignored,  the  overall  pattern  oi  the  remaining  values  is  ;hat  associated  with 
the  semiannual  variation.  The  two  minima  are  in  late  July  1969  and  early 
January  '  '0 ,  and  there  is  a  maximum  between  these  two  date%  at  October* 
November  ,,r>9.  The  v.Lucs  of  density  a*  semi-annual  maximum  and  minima  are 
given  in  (  x'is  2,  The  values  recorded  are  those  sustained  over  10  days 
(excluding  "alues  at  of  geomagnetic  disturbance),  and  are  marked  as 

broken  lines  ir*.  rig.- 

Table  2 


Semi-annual  maximum  and  minima  in  density 
1  l  1969-1970  at  heights  near  240  km 


Central  date 

Values  of 
p245  W8/m3 

July  minimum 

1969  July  25 

0.0520 

October-November  maximum 

1969  October  31 

0.0950 

January  minimum 

1970  January  4 

0.0808 

The  dates  of  the  semi-annual  maximum  and  minima  are  in  agreement  (within  iO  days) 

8 

with  those  obtained  in  the  exosphere  ,  at  heights  near  1000  km,  and  at  heights 
1  2 

below  200  km  *  ,  The  values  in  Table  2  have  not  been  c6rrected  for  solar 

activity:  there  is  no  need  for  correction  of  the  last  two  values,  for  which 
22  — 2 

S,0  7  *  145  *  10  Wm  Hz  ,  but  the  July  minimum  would  be  increased  by 

about  18%  if  corrected  to  the  same  level  of  solar  activity;  with  this  correc- 

3  3 

tion  the  July  value  of  0.0520  ug/m  would  be  increased  to  0.0614  ug/m  . 

5.5  Variation  with  height 

The  values  of  are  plotted  against  y^  in  Fig. 10,  to  give  the 

variation  of  density  with  height.  The  values  plotted  as  circles  are  intended 

to  represent  average  values  of  density  derived  from  the  orbit  of  1969-20B, 

3 

defined  as  those  which  give  values  of  p  between  0.075  ug/m  and 


ifiililiiiiiLiiiiliniiiii. 
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0.100  pg/m  in  Fig. 9.  Values  outside  this  band  are  plotted  as  crosses.  The 

7 

dotted  curve  in  Fig. 10  represents  the  values  given  by  Jacchia  for  exospheric 
temperature  1100  K.  This  is  the  model  tc  which  the  present  values  of  density 
were  scaled  at  heights  between  18G  and  190  km  (see  section  4).  The  broken 
curve  in  Fig. 10  shows  the  Jacchia  model  for  exospheric  temperature  1000  K.  Had 
the  density  values  been  scaled  to  this  curve  (at  the  same  heights),  they  would 
be  reduced  by  8%.  In  the  group  of  values  at  heights  of  220-230  km  an  82 
decrease  would  result  in  more  points  lying  above  the  broken  curve  than  below 
it,  whereas  the  opposite  applies  for  the  dotted  curve:  so  it  appears  that  the 
best  model  would  be  for  exospheric  temperature  of  about  1050  K.  It  can  be 
concluded  that  the  model  fits  the  mean  values  well,  and  that  the  appropriate 
exospheric  temperature,  1050  K,  is  close  to  the  temperature  that  would  be 


expected  to  fit  the  conditions  (S 


DISCUSSION  OF  RESULTS 


10.7 


145). 


The  density  ratios  obtained  from  this  study  for  the  semi-annual  variation 

are  1.55  for  the  October  1969/July  1969  density  ratio  after  the  correction  for 

solar  activity;  and  1.18  for  the  October  1969/January  1970  density  ratio. 

0 

These  values  fit  in  remarkably  well  with  those  obtained  by  Cook  from  the  orbit 
of  1964-63C  at  1070  km,  and  are  plotted  in  Fig. 11,  together  with  values1  from 
1967-31A  at  185  km.  Such  good  agreement  in  the  trend  of  the  variations 
suggests  that  the  factors  influencing  the  density  in  the  exosphere  also  apply 
to  the  atmosphere  around  200  km. 

The  sudden  rise  in  density  about  MJD  40600  apparent  in  Fig. 9  needs 

further  examination,  as  this  is  a  time  when  the  density  should  still  be  fairly 

low  following  the  January  minimum.  Also  the  increase  in  solar  activity  at  the 

time  of  this  large  increase  in  density  is  not  as  great  as  in  several  of  the 

peaks  in  previous  months,  where  there  was  no  such  increase  in  density;  so  the 

increase  in  density  is  rather  surprising.  However  the  analysis  of  the  orbits 

8 

of  1 965-1 6G  and  1964-63C  shows  that  there  was  also  a  large  increase  in  density 
at  this  time  at  heights  near  1000  km.  This  would  seem  tc  confirm  that  the  rise 
in  density  apparent  in  Fig. 9  is  real,  and  not  due  to  a  sudden  change  in  cross- 
sectional  area  of  the  satellite.  This  finding  illustrates  once  again  that  the 
10.7-cm  radiation  is  an  imperfect  index  of  the  extreme  ultraviolet  radiation 
which  influences  the  density,  and  therefore  sometimes  fails  to  shew  the 

correlation  clearly.  This  increase  in  density  at  MJD  40600  was  not  detected 

2 

in  the  analysis  of  1969'-! 08A  because  this  satellite  had  a  lower  perigee 
height  and  the  density  was  evaluated  at  about  160  km,  where  the  atmosphere  is 
much  less  affected  by  solar  activity. 
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CONCLUSION'S 


(1)  The  main  purpose  of  this  paper  is  to  establish  the  variations  in  air 
density  at  heights  of  200-250  km,  between  July  1969  and  February  1970,  Fig. 9 
gives  the  results. 

(2)  It  is  noticeable  that  density  increases  sharply  at  times  of  geomagnetic 
disturbance,  as  expected,  and  that  the  10.7-cm  radiation  energy  is  not  always 
a  good  index  of  the  influence  of  solar  activity. 

(3)  The  July  1969  semi-annual  minimum  is  much  stronger  than  that  in 
January  1970,  and  this  agrees  remarkably  well  with  results  previously 
obtained  in  the  exosphere  at  heights  near  1000  km  (Fig. 10). 
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